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We study the structure of the nuclear force by using a path-integral hadronization 
approach in a chiral quark-diquark model. After the construction of the chiral 
quark-diquark model, we hadronize it to obtain a meson-baryon Lagrangian. The 
effective meson-baryon Lagrangian incorporates chiral symmetry and the com- 
posite description of the mesons and baryons. Using the effective meson-baryon 
Lagrangian we investigate the structure of the nuclear force in the simple case 
of neglecting the axial-vector diquark. It is shown that the meson-baryon La- 
grangian contains two kinds of the nuclear force; the meson-exchange interaction 
and a quark-diquark loop interaction. It is also shown that the quark-diquark loop 
interaction consists of the scalar and vector interactions. The properties of these 
interactions are discussed. 



1 Introduction 

The origin of the short range repulsion in the nucleon-nucleon(iViV) inter- 
action is a long-standing problem. The NN interaction below the meson 
production threshold is phenomenologically well understood, having been fit 
to the phase shift analysis, but its microscopic understanding is still needed. 
While the long range part is well described by the meson exchange picture, 
there are several different approaches to the descrip tion of the short range 
part, including meson exchange and quark exchange I l | ^ l 3 l 4 |5lfil Because in 
the short range region the two nucleons have a substantial spatial overlap, 
it is commonly assumed that the internal structure of the nucleon gives a 
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sizable contribution to the nuclear force at such short distances. To obtain 
a comprehensive understanding of the nuclear force from the short to long 
range distances, models which incorporate both the meson exchange and the 
internal struture of the nucleon are needed. 

In a recent publication, we studied the structure of the nuclear force using 
the path-integral hadronization approach to a chiral quark-diquark model "A 
This method incorporates two aspects of chiral symmetry, which naturally de- 
scribes the pion exchange interaction, and of the internal structure of hadrons. 
This method employs an extended model of the Nambu- Jona-Lasinio type ^ 
with the interactions that are not only of the quark-antiquark type but also 
of the quark-diquark type^l The quarks and diquarks were integrated out to 
generate an effective Lagrangian for mesons and baryons while maintaining 
important symmetries, such as the gauge and chiral symmetries. The hadron 
structure was then described in terms of its constituents: a quark and an 
antiquark for mesons and a quark and a diquark for baryons. It was also 
pointed out that the resulting effective Lagrangian contains various interac- 
tions among hadrons, such as meson-meson, meson-baryon and baryon-baryon 
interactions. 

In our model, all components of the NN force are contained in an effec- 
tive Lagrangian that is written in a concise form as a trace-log. Then, the 
expansion of the trace-log terms produces an NN force that is described as 
meson exchanges at long and medium ranges and quark-diquark exchanges at 
short ranges. The latter was then shown to contain various types of non-local 
interactions, including a scalar iso-scalar type and a vector iso-scalar type. 
We evaluated the properties of the scalar and vector type interactions; their 
ranges, effective masses and strengths. 

In this paper, we report the results of the study of the nuclear force in the 
simple framework where the axial-vector diquark is neglected. In the previous 
work, we suggested that the neglecting the axial- vector diquark does not affect 
to the size of the nucleon, hence we can properly evaluate the ranges of the 
NN interactions with only scalar diquarks. 

We organize this paper as follows. In § 2, we briefly give the derivation 
of the trace-log formula in the path-integral hadronization of the NJL model 
with quark-diquark correlations. In § 3, terms containing the NN interaction 
are investigated in detail, where the general structure of the NN amplitude 
is presented. We present a sample numerical calculation for the case in which 
there is only a scalar diquark. The present study of the NN interaction is not 
quantitatively complete, but it will be useful in demonstrating some important 
aspects of the nuclear force, in particular that the range of the short range 
interaction is related to the intrinsic size of the nucleon. The final section is 
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devoted to a summary. 



2 Effective Lagrangian for mesons and nucleons 

We now briefly review the method to derive an effective Lagrangian for mesons 
and nucleons from a quark and diquark model of chiral symmetry following 
the previous work of Abu-Raddad et alii. We start from the NJL Lagrangian, 

G 

%JL = <?W ~ m o)<l + 2 [^ g ) 2 + (^75T<?) 2 ] . (1) 

Here, q is the current quark field, f represents the isospin (flavor) Pauli matri- 
ces, G is a dimensional coupling constant, and mo is the current quark mass. 
In this paper, we set too = i.e. we work in the chiral limit. As usual, the 
NJL Lagrangian is bosonized by intro ducing me son fields as collective auxil- 



iary fields in the path- integral method* * * \ At an intermediate step, we 
find the following Lagrangian: 

C qaw = q (ip - (a + i75 f • tt)) q - j^{a 2 + 7? 2 ) . (2) 

Here a and 7? are scalar-isoscalar sigma and pseudoscalar-isovector pion fields, 
as generated from a ~ qq and 7? ~ iqr^^q, respectively. For our purp ose, it 
is convenient to work in a non- linear basis rather than the linear one ^3Jl2 1 
This is realized through the chiral rotation from the current (q) to constituent 
(x) quark fields: 

X = tsq, & = I J ) , (3) 

where f 2 = a 2 + tt 2 . Thus, we find 

c' x ^ = x ^-f-^-Mx-^f\ (4) 

where 

= -\ + dptf) , a, = - l - (d^d ~ d„tf) (5) 

are the vector and axial-vector currents written in terms of the chiral field, 

/ I * — * — *\ 1/2 

t -(=±p) ■ (6> 

The Lagrangian describes not only the kinetic term of the quark, but also 
quark-meson interactions such as the Yukawa one, Weinberg- Tomozawa one 
etc. 
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In the model we consider here, we introduce diquarks and their interaction 
terms with quarks. We assume local interactions between quark-diquark pairs 
to generate the nucleon field. As suggested by a method of constructing a 
local nucleon field, it is sufficient to consider two types of diquarks, a scalar, 
isoscalar diquark, D, and an axial-vector, isovector diquark, In this 

work, we consider only the scalar diquark. Hence, our microscopic Lagrangian 
for quarks, diquarks and mesons is given by 

C = X(ip ~f-i>- H)x - 7^f + D\d 2 + M 2 S )D + G X ^D X ■ (7) 

In the last term, G is a coupling constant for the quark-diquark interaction. 

Now, the hadronization procedure can be carried out straightforwardly by 
introducing the baryon fields as auxiliary fields, B ~ D\, and by eliminating 
the quark and diquark fields in Q. The final result is written in a compact 
form as 

£cff = -^/ 2 - i tr ]n(iP -J-i>- ^75) - 4 BB + i tr ln(l - S) . (8) 

Here, the trace is taken over space-time, color, flavor and Lorentz indices, and 
the operator S is defined by 

S = A T B SB, (9) 

In this equation, S = (if) — f — ji — ^7 5 ) _1 , A = (d 2 — AfJ) -1 are the propa- 
gators of the quark and scalar diquark, respectively, and transposed diquark 
propagators, as denoted by the superscript T, are employed. Note that the 
quark propagator S contains the interactions with pions. Through this in- 
teractions, the nucleon-pion interactions are obtained. Though the effective 
meson-nucleon Lagrangian (JHJ looks simple, it contains many important phys- 
ical ingredients when the trace-log terms are expanded: 

• It generates a meson Lagrangian in a chirally symmetric manner. Up to 
fourth order in the meson fields, it produces precisely the Lagrangian of 
the linear sigma model with the realization of the spontaneous breaking 
of chiral symmetry. Hence, the vacuum expectation value / turns out to 
be the pion decay constant f v . 



From the second trace-log term, a nucleon effective Lagrangian is derived. 
In a previous paper, the kinetic term of the nucleon was investigated, and 
the mass of the nucleon was computed at the one- loop level 



In the nucleon effective Lagrangian, meson-nucleon couplings appear 
through the diagrams, as shown in Fig. ^ Their strengths and form 
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Figure 1. A quark-diquark one-loop diagram for the meson- nucleon 
Yukawa vertex. The solid, double, dotted and triple lines represent 
the quark, diquark, meson and nucleon, respectively. The blobs repre- 
sent the three point quark-diquark-baryon interaction. 




Figure 2. A loop diagram for the NN interaction (left) and the equiv- 
alent diquark exchange diagram. 



factors can be computed with the use of the underlying quark-diquark 
dynamics. Using these vertices, meson-exchange interactions are con- 
structed, the photon-nucleon couplings are also obtained in the same 
manner. The radii and magnetic moments were evaluated without the 
axial- vector diquarlPd. 



• There are diagrams that contain many nucleon fields. For instance, NN 
interactions are expressed as one- loop diagrams, as shown in Fig. EI This 
term describes the short range part of the NN interaction. In this paper, 
we focus our attention mostly on the NN interaction derived from the 
one-loop diagrams. 
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3 Short range interaction 



We evaluate the short range interaction described by the quark-diquark loop, 
as shown in Fig. [3 Using the interaction vertices given in Eq. ©, it is 
straightforward to compute the amplitude for the quark-diquark loop: 



d 4 k 

M NN = ~iN c Z 2 



(2tt)4 

Bjp'^ + m q )B( Pl )B(p' 2 )^ 2 -A + I + m q )B{p 2 ) 
[(pi - fc) 2 - MJP 2 - mm P [ - fc) 2 ~ Af|][(p 2 + fc) 2 - m 2 ] 



.(10) 



where N c is the number of colors and Z is the wave-function renormalization 
constant. In this equation, the momentum variables are assigned such that p\ 
and p[ (p2 and p' 2 ) are for the pair of contracted baryon fields, B(p[) ■ ■ ■ B(p\) 
(B(p' 2 ) ■ ■ ■ B(jp2)), and the momentum transfer is defined by q — p' x — p\. 
Note that the momentum q is carried by the diquark pair. The amplitude 
defined in this way can be interpreted as a direct term in the local potential 
approximation. When computing physical quantities such as phase shifts, 
we need to include the exchange term that is obtained by interchanging the 
momentum variables p[ «-> p' 2 . Although the one-loop integral (|1(J|I converges 
when the scalar diquark is included, we keep the counter terms of the Pauli- 
Villars regularization. Because our model is a cut-off theory with a relatively 
small cutoff mass, Apy = 0.63 GeV, the counter-terms play a significant role. 
However, because we include only the scalar diquark, we do not attempt to 
make a comparison at the quantitative level. Rather, in the following, we 
study some basic properties of the amplitude itself, mostly the interaction 
ranges extracted from Eq. l|lUfl . 

Let us evaluate the integral in the center-of-mass system for elastic scat- 
tering: 

Pi = (Ep, p), P2 = (Ep, -p), 

p[ = (E p -,, p'), p' 2 = [Ep,, -p') , |pl = |p"| . 

To proceed, we write the amplitude (|10|l as 

M NN = F S (P, q)(BB) 2 + F V (P, q)(B^B) 2 + • • • , (11) 

where 

P = p'+p, q = p'-p. (12) 

Roughly speaking, the g-dependence expresses the interaction range in the 
i-channel, whose Fourier transform is interpreted as the r-dependence of the 
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local potential, while the P-dependence expresses the non-locality of the in- 
teraction. 

In Eq. , we have defined the coefficients F$ and Fy as scalar and vector 
interactions, respectively. The ellipsis in Eq. (|II|I then include terms involving 
external momenta pi, such as (BT(pi)B)(BB) and (BB)(BT(pi)B), where 
r(pi) are 4x4 matrices involving pi. These momentum dependent terms are, 
however, expected to play a less important role than the dominant components 
of the scalar and vector terms, as in boson exchange models. Therefore, we 
consider here the corresponding terms of scalar and vector types. It turns out 
that the interaction coefficients defined in this way are attractive for Fs and 
repulsive for Fy, 

As anticipated, the amplitude is highly non-local, as the quark-diquark 
loop diagram implies. As a matter of fact, the interaction range for the diquark 
exchange in the t-channel is shorter than that for the quark exchange in the 
s-channel. This is shown explicitly in Fig.El Nevertheless, we proceed further 
and carry out an expansion in P. We obtain 

F i (P,q)=F i (P,q)\ p ^ +P~F i (P,q)\ p=0 + ... , (13) 

oP 

where i stands for S or V. The resulting q dependent function, in particular 
the first term, can be interpreted as the Fourier transform of a local potential 
as a function of the relative coordinates r = x\ — xi. Then, the general struc- 
ture of the NN potential can be studied by performing the non-relativistic 
reduction of the amplitude, Eqs. (|ll|l and ljl3JI . It contains central, spin-orbit 
and tensor components that accompany functions of non-locality P. 

We now discuss the functions Fs and Fy in Eq. i|ll|) for the leading-order 
terms of Eq. (|13fl . Because we cannot write the resulting q dependent functions 
in an analytic form, we have carried out numerical calculations employing 
several different values of G parameter for different binding energies and the 
size of the quark-diquark bound state. 

At this point, we explain our parameters and the regularization scheme. 
We follow the scheme presented in Ref. ^ except for the treatment of G. 
The A py is the cutoff mass of the Pauli-Villars method used in this work to 
regularize divergent integrals. We note that our model is unrcnormalizable. 
The NJL coupling constant G and the cut-off mass A py are fixed to generate 
the constituent quark mass m q an d the pion decay constant through the NJL 
gap equation in the meson sector The mass of the scalar diquark 

Ms is determined in t he NJL model by solving the Bethe-Salpeter equation 
in the diquark channel I Then, we have one free parameter, the quark- 
diquark coupling constant G. The strength of G controls the binding or size 
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q- or P 2 [GeV 2 ] q 2 or P 2 [ Ge V 2 ] 

Figure 3. Normalized scalar and vector form factors as functions of q 2 = 
\q 2 \ with fixed P = (solid curve) and those of P 2 — \P 2 \ with 
fixed q — (dashed curve). The calculations were performed using the 
parameter set given in Table 1. 



Table 1. Model parameters. 



Mjv [GeV] 


m, [GeV] 


M s [GeV] 


A PV [GeV] 


G [GeV" 1 ] 


( r 2}l/2[ fm ] 


set I 0.94 


0.39 


0.60 


0.63 


271.0 


0.77 


set II 0.85 


0.39 


0.60 


0.63 


445.9 


0.54 



of the nuclcon and generates the mass of the nucleon M^. Then we perform 
numerical calculations by varying the coupling constant G. In Tabled we list 
typical parameter values as used in Ref. 1-4 

First, we determine the strengths of the interactions by extracting cou- 
pling constant squares gf as 

1^(0,0)1 = -^, 5*>0, (14) 

where the effective meson masses rm are evaluated using the inverse of (I15|l 
and are plotted in Fig. El The results are displayed in Fig. 0] as functions of 
the size of the nucleon, (r 2 ) 1 / 2 . The coupling strengths can be compared with 
the empirical values gs ~ 10 and gy ~ 13 El The present results are strongly 
dependent on (r 2 ). When only a scalar diquark is included, the scalar interac- 
tion becomes much stronger than the vector interaction. Phenomenologically, 
the vector (omega meson) coupling is stronger than the scalar (sigma meson) 
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coupling. 

We should comment on the effect of the axial-vector diquark. The loop 
integral containing the axial- vector diquarks diverges and has a large numer- 
ical value even when it is regularized. Therefore, it significantly affects the 
absolute values of the loop integrals as well as that of the baryon self energy 
which is necessary to extract the normalization factor Z . Therefore, it is ex- 
pected that the strengths are affected significantly, but the ranges and the 
corresponding masses, which are computed from the momentum dependence 
of the normalized loop integrals, are not affected much. 

80.0 
60.0 
% 40.0 
20.0 



""0.4 0.6 0.8 1.0 1.2 

</-V n (fin) 

Figure 4. Scalar and vector coupling constants as functions of the nu- 
cleon size, (r 2 } 1//2 . 

Let us discuss the interaction ranges. In Fig.|31 we show the g-dependence 
of the zeroth order coefficients in Eq. for three different sizes of the 

nucleon. It is obvious from Fig. [S] that as the size of the nucleon becomes 
smaller the interaction ranges become shorter. More quantitatively, we define 
the interaction ranges Ri by 

(15) 

q 2 ^0 

which are related to the mass parameters of the interaction ranges as m, = 
y/E/Ri. It is interesting that the range (and hence mass) parameters of the 
interactions are approximately proportional to the size of the nucleon (r 2 ) 1 / 2 , 
as shown in Fig. When the parameter set I listed in Table ^ is used, the 
mass parameters are about 650 MeV and 800 MeV for the scalar and vector 
interactions, respectively, which are very close to the masses of the sigma and 
omega mesons. If, however, we use the parameter set II for a nucleon size of 
about 0.5 fm, then the two masses become ms ~ 800 MeV and my ~ 1000 




R? = -6 



1 dF t 
FM 2 ) dq 2 
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MeV. 




0.0 0.5 1.0 0.0 0.5 1.0 

q 2 lGeV 2 ] cf [GeV 2 ] 

Figure 5. The scalar (left panel) and vector (right panel) form factors 
Fs,v (q 2 , P 2 = 0)/F s ,v{q 2 = 0, P 2 = 0), (form factors are normalized 
to 1 at q 2 = 0). The curves correspond to A'/jv=0.98, 0.85, 0.70 GeV 
from bottom to top. 

In the present analysis, apart from the absolute values of the interaction 
strengths, the interaction ranges for the scalar and vector interactions have 
been produced appropriately with the quark-diquark loop diagram. One ques- 
tion concerns the small but non-negligible difference between the ranges in the 
scalar and vector channels, which is consistent with empirical results. This 
can be roughly understood from the dimensionality of the loop integral. As 




0.4 0.6 0.8 1.0 1.2 0.4 0.6 0.8 1.0 1.2 

<r > (fm) <r > " (fm) 

Figure 6. The interaction range Ri (left panel) and the corresponding 
mass parameter nn (right panel) as functions of the nucleon size (dis- 
tribution of a quark and a diquark). 
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seen from Eq. 1|1U[) . the integrand for the vector interaction is of higher order 
with respect to the loop momentum than the scalar interaction. Because of 
this, the vector part reflects shorter distant dynamics and produces a shorter 
interaction range. 

Although we have obtained reasonable results for the interaction ranges, 
the short range part of the nuclear force must be repulsive. Hence, current 
approach is not satisfactory. We believe that the inclusion of the axial- vector 
diquark will give a sizable contribution to the short range repulsion because 
the quark-axial- vector diquark-baryon vertex is basically spin-spin type inter- 
actions, and thus correct this shortcoming of our present approximation. 

4 Summary 

In this paper we have studied the NN interaction using a microscopic theory 
of quarks and diquarks. Nucleons were described as quark-diquark bound 
states. The quark and diquark degrees of freedom were integrated out using 
the path-integral method, and an effective Lagrangian was derived for mesons 
and baryons. The resulting trace-log formula contains various meson and 
nucleon interaction terms, including the NN interaction in the short range 
region expressed as a quark-diquark loop. Hence, the NN interaction can be 
naturally expressed as meson exchanges at long ranges and quark or diquark 
exchanges at short ranges. 

For short range interactions, we have computed a quark-diquark loop cor- 
responding to the direct term of the two-nucleon interaction. The resulting 
interaction is highly nonlocal. We have extracted scalar and vector type in- 
teractions in the local potential approximation. It turns out that the scalar 
term is attractive, while the vector is term repulsive. 

In the present paper our numerical calculations contained only scalar 
diquarks, as a first step toward a full calculation. Therefore, we have con- 
centrated our study mostly on the interaction ranges or, equivalently, the 
masses, because the magnitude of strengths is affected by the axial- vector di- 
quark. On the other hand, the interaction ranges were better studied in the 
present study including only scalar diquarks, reflecting the size of the nucleon. 
Consequently, the mass parameters of the interaction ranges for the scalar and 
vector interactions were found to be about 650 MeV for the scalar type and 
about 800 MeV for the vector type, once again with the nucleon size set at 
0.77 fm. Hence in our model, the scalar-isoscalar interaction emerges from the 
quark-diquark loop at an energy scale similar to that of the sigma and omega 
meson exchanges. The existence of the two components in the nuclear force, 
boson exchanges and quark-diquark exchanges, is a general feature when we 
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consider a model of nucleons that is composed of a quark (and diquark) core 
surrounded by meson cloud. 

The present result encourages us to further study baryon properties by 
extending the model to include the axial-vector diquark. We have already 
started such a study. When the axial- vector diquark is included, loop integrals 
diverge more strongly than in the case of the scalar diquark, due to the massive 
vector nature of the propagator. Although this causes the numerical study to 
be more complicated than in the present case, this work is in progress. 
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